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ABSTRACT 

In this paper we obtain constraints on the cosmological density parameters, 
Q m and fib, by comparing the preliminary measurement of the QSO power 
spectrum from the two degree field QSO redshift survey with results from an 
analytic technique of power spectrum estimation, described in this paper. We 
demonstrate the validity of the analytic approach by comparing the results 
with the power spectrum of an N-body simulation. We find a better fit to the 
shape of the QSO power spectrum for low density models with fi m = 0.1 — 0.4. 
We show that a finite baryon fraction Qt,/Q, m = 0.2, consistent with obser- 
vations of the CMB anisotropies and nucleosynthesis, fits the observational 
result of the 2QZ survey better, though the constraint is not very tight. By 
using the Fisher matrix technique, we investigate just how a survey would 
be required before a significant constraint on the density parameters can be 
made. We demonstrate that the constraint could be significantly improved if 
the survey was four times larger. 
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1 INTRODUCTION 

Since the pioneering work of Osmer in the 1980's (Osmer 1981), QSOs have been used to 
probe the high redshift Universe. It was soon shown that QSOs were clustered (Shaver 1984; 
Shanks et al. 1987; Iovino & Shaver 1988) but results were limited due to the relatively small 
QSO samples available, see Croom & Shanks (1996) and references therein for a summary 
of early QSO clustering results. The two degree Field QSO Redshift Survey (2QZ) will 
dramatically improve things as it will contain at least a factor of 25 more QSOs than than 
previous QSO surveys. The 2QZ group has recently reported preliminary clustering results 
based on an initial sample of 10,000 QSOs (Croom et al, 2001, Hoyle et al. 2002; 2001, 
Outram et al. 2001). 

Observational results and theoretical predictions are most directly compared through the 
use of numerical simulations. Indeed, the 2QZ group has utilized one of the Hubble Volume 
simulations run by the Virgo consortium (Frenk et al. 2000) to interpret their results. It is a 
huge N-body simulation of horizon box size, containing 1 billion mass particles output along 
a light-cone. However, numerical simulations of this size are slow and expensive to run. A 
simple, semi-analytic formula that reproduces numerical results but which doesn't require 
intensive computation would be useful as it allows a wide range of parameter space to be 
explored in a short time frame. Fortunately, several examples appear in the literature (e.g. 
Yamamoto et al. 1999; Suto et al. 2000; Hamana et al. 2001b; and references therein). 

In this paper, we describe such an analytic approach. We demonstrate the validity of 
this technique by comparing the semi-analytic formula with results from N-body simulations 
and mock 2QZ samples. As a demonstration of the usefulness of this approach, we use the 
analytic formula to place constraints on the cosmological density parameters by comparing 
the theoretical predictions to the 2QZ power spectrum of Hoyle et al. (2002). By using the 
Fisher matrix approach we also make predictions for the size of the survey that would be 
required before the density parameters can be tightly constrained, (see e.g., Tegmark 1997, 
Tegmark et al. 1998) 

This paper is organized as follows: In section 2, we briefly summarize a theoretical formula 
for the power spectrum that incorporates the light-cone effect, the linear redshift-space 
distortions and the geometric distortion. In section 3, we compare the analytical power 
spectrum formula with power spectra results from the N-body simulation and mock QSO 
samples. In section 4 we place constraints on the density parameters, fl m and Qb, by fitting 
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the analytic formula to the QSO power spectrum measured from the 10k 2QZ sample. In 
section 5 we make predictions for how well the density parameters may be constrained with 
yet larger surveys and Section 6 is devoted to summary and conclusions. Throughout this 
paper we use the unit in which the light velocity c equals 1. 



2 ANALYTIC FORMULAS 

In this section we briefly review an analytic formula for the spatial power spectrum of 
cosmological objects. The analytic form of the power spectrum has to include light-cone 
effects, redshift-space distortions and geometric distortions. We follow, for example, Suto 
et al. (2000) and define 



dzW(z)P a (k,z) 

P hC {k) = 1 t , (1) 

dzW(z) 



where 

dN/dz is the number density of objects per unit redshift per unit solid angle, s = s(z) is the 
distance-redshift relation defined by equation (||), and P£ (k,z) is the local power spectrum 
defined on a constant time hypersurface at redshift z. Thus the light-cone effect is taken 
into account by averaging the local power spectrum Pq(/c, z) over the redshift range. The 
distance s = s(z) is the comoving distance which, in a specially flat universe, is given by 
1 f z dz' 

s(z) = TT / / (3) 



where Ho = 100/ikm/s/Mpc, the Hubble parameter, and Q m is the total matter density 
parameter. In the present paper we fix Q m = 0.3 and h = 0.7 in equation (|3|), consistent 
with the values used in the analysis by Hoyle et al. (2002). For the local power spectrum, 
P a (/c, z), we model 

P Q (k, Z) = — — ~. t^t 

c^zjc^z) 2 

x / dfiP QS o(q\\ ~^—,q±^ k ^ 1 ~^ ,z), (4) 

where Pqso(<1\\ iQ±, z ) is the QSO power spectrum, q ]{ and q ± are the wave number components 
parallel and perpendicular to the line-of-sight direction in the real space, and we define 
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where r(z) is the comoving distance of the universe, found using equation (3) but with 
varying Q m . In the present paper, we only consider the spatially flat cosmological model 
with a cosmological constant. Because an analysis of the observational data is performed by 
fixing the cosmological model, i.e., by equation (3), we need to take the geometric distortion 
effect into account when comparing the observational result with models using different 
cosmological parameters. The geometric distortion effect is included by scaling the wave 
number in equation (4) by the factors C\\(z) and c ± (z). 

We model the QSO power spectrum of the distribution with the linear distortion (Kaiser 
1987) by 

P Q So(q }l ,q±,z) = (l + ^-i) 2 b(z) 2 P m&ss (q,z), (7) 



v b(z) q 2 / 

where q 2 = q 2 + q 2 ±) b(z) is a scale independent bias factor, P maBB (q, z ) is the CDM mass 
power spectrum, and we defined f(z) = d In D (z)/d In a(z) with the linear growth rate D(z) 
and the scale factor a(z). For a phenomenological correction for the nonlinear velocity field, 
a damping function might be multiplied on the right hand side of equation fl7|). This is a 
very small correction on the scales we consider in this work so we neglect it here, (see e.g., 
Mo, Jing & Borner 1997; Magira, Jing & Suto 2000) 

The two-point correlation function, corresponding to the spectrum (1), can be written 

e LC (i?) = — / dkk 2 3o (kR)P Lc (k). (8) 

2,71 Jo 

Formula ([I]) is obtained under the short distance approximation, R(= 2n/k) <C L, where L 
is the size of a survey area. Hence the validity of the use is limited to small length scales. 
The finite size effect means we cannot obtain a robust estimation of the correlation function 
or the power spectrum on length scales larger than the size of a survey. For the 2QZ survey, 
the limit is L ~ 10 3 /i -1 Mpc. For the correlation function (§), it was reported that the 
analytic expression is in good agreement with results of the N-body numerical simulation 
(Hamana et al. 2001a; 2001b), though the comparison might be limited to a finite range of 
length scales. 

We only consider the linear power spectrum as nonlinear effects are negligible on the 
length scales on which the power spectrum from QSOs can be measured. Hence we adopt 

P ma ss(g, z) = Aq n T(q, Q m , Q b , h) 2 D(z) 2 , (9) 
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where A is a normalization constant, T is the transfer function, and n is the index of the 
initial spectrum, for which we assume n = 1 in the present paper. We adopt the fitting 
formula of the cold dark matter transfer function by Eisenstein & Hu (1998) which is robust 
even when the baryon fraction is large. 

3 COMPARISON WITH N-BODY NUMERICAL SIMULATION 

In this section we investigate how well the analytic power spectrum reproduces a power 
spectrum measured from an N-body numerical simulation and a power spectrum measured 
from mock 2QZ catalogues. 

The simulation is one of the Hubble Volume (HV) simulations. It has a ACDM cosmology 
with fi m =0.3, £7^=0.04 and ^a=0.7, h—0.7 and normalized to o"8=0.9 at present day and 
with a shape parameter of T = 0.17. The simulation has a light-cone output such that the 
evolution of the dark matter is fully accounted for. It extends to redshift ~ 4, and covers an 
area of 15x75deg 2 , which we split into three 5x75deg 2 strips to match the geometry of the 
2QZ strips. 

For the first comparison, we simply take the dark matter points in real space, over the 
range 0.3 < z < 2.5. We apply the selection function of the 2QZ survey shown, for example, 
in Hoyle et al. (2002), to the points and then measure the power spectrum. We do this in each 
of the three 5x75deg 2 strips and average the result together. The analytic power spectrum 
assumes the same dN/ dz distribution, that the points are in real space and assumes the same 
cosmology. Figure 1 shows a comparison of these two power spectra. The solid and dotted 
lines, respectively, show the results of the analytic prescription and the numerical simulation. 
There is good agreement (within ~ 10 %) on scales COl/iMpc" 1 < k < 0.1/iMpc" 1 . 

As a second comparison, we compare the analytic power spectrum to the power spectrum 
of mock QSO catalogues drawn from the HV. The mock catalogues are constructed as follows: 
First, the selection function of the 2QZ is applied to the dark matter particles in redshift- 
space over the redshift range 0.3 < z < 2.5. Then the dark matter particles are biased such 
that the clustering approximately matches that of the 2QZ. It is assumed that the QSO 
clustering amplitude remains constant with redshift, consistent with the results of Croom et 
al. (2001), when the mock catalogues are constructed. The points are then sparse sampled 
to match the number density of the 2QZ. Again the power spectrum is measure from each 
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strip and the results are averaged together. Full details on the construction of the mock 
catalogues are given in Hoyle (2000). 

For this comparison, the amplitude of the analytic power spectrum has to be determined. 
To do this, we minimizing x 2 defined by equation fllPp. 

*= V aW.P (10) 

where P sim (ki) is the value of the power spectrum from the simulation at hi, AP(ki) 2 is 
the variance of errors in Figure 2, and N is the number of the data points. In the analytic 
modeling of the clustering bias, we assumed the form 

Hz) = T ^- r (11) 
v ; [D(z)]p v ; 

where bo and p are the parameters. We determine bo to minimize x 2 f° r each pair of Q m and 
fib. Croom et al (2001) report that the QSO clustering amplitude does not show significant 
time-evolution, therefore we set p = 1. However, this choice does not alter our conclusions. 
For example, when we set p = 2 the values of x 2 change by only a few percent. This confirms 
the result that the power spectrum is only sensitive to a mean amplitude of the bias and 
is insensitive to the speed of the redshift evolution (Yamamoto & Nishioka 2001). Figure 2 
shows that a more realistic numerical result can still be fit by the analytic approach even 
though the bias is calculated in a different way. 

Figure 3 shows the x 2 contours found when different values of the cosmological density 
parameters Q m — Qb/Q m are used in the analytic formula. Each panel corresponds to the 
constraint using (a) the 17 points (N = 17), (b) the 16 points, (c) the 15 points and (d) the 
14 points, respectively, of the data from left to right in Figure 2. The cross point in each 
panel is the target parameter of the simulation. In the case where we including data points 
of small length scales, the minimum of x 2 is slightly inconsistent with the target parameter. 
The disagreement on these small scales k ~ 0.2 can not be explained by including the 
nonlinear corrections of the density perturbations and the finger of God effect. It might be 
inferred that long mean separations between QSOs cause this feature. In general, however, 
agreement of the numerical simulation and the analytic approach is acceptable over the 
range of scales 0.01/iMpc" 1 < k < 0.1/iMpc -1 . 
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4 COMPARISON WITH AN OBSERVATIONAL RESULT 

As an application, we compare the power spectrum measured from the 10k 2QZ dataset to 
the analytic formula as a way to constrain cosmological parameters. We simply define \ 2 m 



data, for which we use the result of Hoyle et al. (2002). Figure 4 shows the observational 
data points from Hoyle et al. (2002) compared to the best fit analytic power spectrum. 
Figure 5 displays the contours of \ 2 f° r various cosmological models on the fl m — Q b /fl m 
plane. Each panel shows the contour using the data of (a) the 17 points, (b) the 16 points, 
(c) the 15 points and (d) the 14 points, respectively, from left to right in Figure 4, where 
the wave numbers of the observational data used in Figure 5 are same as those in Figure 
3. From Figure 5 it is clear that the QSO power spectrum favors a low density universe 
with Q m = 0.1 — 0.4 rather than the standard CDM model with Q m = 1. The minimum 
of the x 2 is located at Q m ~ 0.2 and Qb/Q m ~ 0.2. The dashed curve shows = 0.04, 
inferred from results of the comic microwave background (CMB) anisotropics and the big 
bang nucleosynthesis (BBN) (see e.g., Turner 2001). Our result Q m ~ 0.2 is consistent with 
the previous result of the 2QZ power spectrum, which reports the best fit value Vt m h = 0.1 
(Hoyle et al. 2002, cf. Outram et al. 2001), but, is smaller than Q m ~ 0.3, which is obtained 
from the large scale structures of galaxies (e.g., Peacock et al. 2001). It is unclear whether 
the discrepancy has any physical meaning or not, because the 1 — a contour is broad. An 
interesting fact is that the QSO power spectrum is better explained with the finite baryonic 
component, though the peak of \ 2 is broad and thus the constraint is not that tight (cf. 
Miller et al. 2001, Peacock et al. 2001). 

5 TIGHTNESS OF THE CONSTRAINT 

It is useful to investigate how much larger the 2QZ Survey would have to be before a tight 
constraint on the density parameters can be obtained. For this purpose, we adopt the Fisher 
matrix approach. The Fisher matrix analysis provides a technique by which statistical errors 
on parameters from a given data set can be estimated. In general, the Fisher matrix is defined 



the similar way as equation ( |I~0D but replace P mra {ki) with P' 



(ki), i.e., the observational 



by 




(12) 
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where L is the probability distribution function of a data set, given model parameters. 
Here, we follow the work of Tegmark (1997). Using the Fisher matrix approach, Tegmark 
investigated the accuracy with which cosmological parameters can be measured from galaxy 
surveys. Following Tegmark (1997), we can write the Fisher matrix element, 



where Af2 is the solid angle of the survey area and 0, denotes the theoretical parameters. 
In deriving equation (|13|), we use the fact that nPqso{k) ~ O(0. 1) where n is the comoving 
mean number density of the QSOs. Here we set fc max = 0.2/iMpc -1 and fc m ; n = O.Ol/iMpc -1 . 

By using the Bayse theorem, the probability distribution in the parameter space can be 
written P(8i) oc exp[— YUj{^i~ 6t r )^ij(@j ~~@f)/2}, where 9j r is the target model parameters. 
Here, the Fisher matrix is the 2x2 matrix, Q m and Qf, matrix. We show how accurately 
these two parameters can be determined for the final 2QZ sample and for a sample that 
is four times larger in Figure 6. Curves in each panel indicate the 68%, 95%, and 99.7% 
confidence region on the Q m — f2&/f2 m plane. Both panels assume the same dN/dz as the 
10k 2QZ catalogue but the left panel assumes Af2 = 750 deg. 2 and the right panel assumes 
AQ = 3000 deg. 2 . This fi gure shows an example of the tightness of the constraint possible 
from the QSO spatial power spectrum alone. Note that the constraint will be tighter if we 
combine the results with other constraints, such as those from the Alcock-Paczynski test 
(Alcock & Paczynski 1979; Outram et al. 2001; Hoyle et al. 2002). 

6 SUMMARY AND CONCLUSIONS 

In this paper we have shown that an analytic approach can accurately reproduce the power 
spectrum of an N-body simulation and mock QSO catalogues over the range of length scales 
k ^ 0.1/i -1 Mpc. By performing a simple x 2 test we have also shown that the observational 
QSO power spectrum is consistent with a simply biased mass power spectrum based on the 
popular CDM cosmology with a cosmological constant. We find that the analytic formula 
better matches the 2QZ power spectrum if the baryon fraction is around 20% of the mass 
fraction. The constraint is not very tight because our analysis is based on just the initial 
sample of 10,000 QSOs from the 2QZ survey. Accumulation of a number of samples will 




(13) 



where we define 




(14) 
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improve the statistical errors. Tighter constraints will be possible if we compare the results 
with other cosmological constraints such as those obtained from the A-test. 
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Figure 1. Comparison of dark matter power spectra in real space with the 2QZ selection function applied. The dotted curve 
shows the result from the numerical simulation and the solid curve shows the result from the analytic approach. The dashed 
curves are the power spectra P a (k, 2 = 0) (the upper curve) and P a (k, z = 2.5) (the lower curve). 
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Figure 2. Comparison of the power spectra of biased particles in redshift space with the 2QZ selection function applied. The 
points with error bars show the mock catalogue power spectrum from Hoyle et al. (2002) and the solid curve shows the result 
using the analytic approach. 
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Figure 3. The contours of \ 2 show the level to which the power spectrum of the numerical simulation and the analytic 
formula agree. The number of points used in calculating \ 2 is shown in each panel. These are the N points from the left shown 
in figure 2. The cross point is the target parameter of the simulation. 
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Figure 4. The observational data, from Hoyle et al. (2002), which we used in constraining the cosmological parameters. The 
solid curve shows the analytic formula using the best fit parameters f2 m = 0.2 and Cl b = 0.04. 
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Figure 5. Contours of \ 2 found by comparing the observed QSO power spectrum with the analytic formula. The format of 
this figure is the same as Figure 3. The number of data points, N, and the wave numbers of the observed power spectrum, 
ki, are the same as those in Figure 3. The dashed curve shows Q(, = 0.04, which is consistent with results from the cosmic 
microwave background anisotropies and big bang nucleosynthesis, (see e.g., Turner 2001). 
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Figure 6. Confidence regions from data sets of QSO samples. The curves show the 68%, 95% and 99.7% confidence regions. 
The left panel assumes a solid angle of Afl = 750 deg. 2 and the right panel assumes a solid angle of Af! = 3000 deg. 2 for the 
survey region. In both cases, we assume the dN/dz of the 2QZ survey, and we set fc max = 0.2/iMpc -1 and fc m i n = O.OlAMpc - 1 . 
The target parameters are f2 m = 0.2 and Hj, = 0.04. 
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